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ABSTRACT 

We study relativistic self-dual Chern-Simons-Higgs systems in the presence of uniform back- 
ground fields that explicitly break CTP. A rich, but discrete vacuum structure is found when 
the gauge symmetry is spontaneously broken, while the symmetric phase can have an infinite 
vacuum degeneracy at tree level. The latter is due to the proliferation of neutral solitonic states 
that cost zero energy. Various novel self-dual solitons, such as these, are found in both the 
symmetric and the asymmetric phases. Also by considering a similar system on a two-sphere 
and the subsequent large sphere limit, we isolate sensible and finite expressions for the conserved 
angular and linear momenta, which satisfy anomalous commutation relations. We conclude with 
a few remarks on unresolved issues. 
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1 Introduction 



In Abelian Maxwell-Higgs systems, there are well-known topologically stable vortex solutions in 
three dimensions. With a specific coupling constants, the energy functional can be saturated by the 
so-called self-dual configurations satisfying certain first-order differential equations When the 
kinetic term for the gauge field includes Chern-Simons term, the corresponding self-dual models 
are also found and studied [Q, ^, ^] . 

An interesting generalization of such self-dual models arises when the system is coupled to 
an external background charge density or an external magnetic field. In particular, Maxwell- 
Higgs systems with the uniform external electric charge density was argued to describe the real 
superconductor more closely than ones without. The self-dual limit of these systems has been 
investigated extensively by one of the authors ||5|. In this paper, we want to consider two more self- 
dual models, which incorporate a Chern-Simons term. A nonrelativistic version with an external 
magnetic field was studied before as an effective field theory for the fractional quantum Hall effects 

i- 

There are several novel features in (Maxwell) Chern-Simons-Higgs systems with the background 
charge density. Their structure is much richer than the systems without. First, homogeneous 
ground states are possible in two rather different manners: either a symmetric phase with a uniform 
magnetic field or an asymmetric phase with a uniform Higgs charge density. (But they obviously 
belong to two different superselection sectors, even if the spatial volume is finite.) Second, a 
magnetoroton mode is possible in the asymmetric phase for certain parameter ranges. This mode 
has the lowest energy at a nonzero wavelength. Third, such a roton mode, when its energy is 
imaginary, leads to the instability of the homogeneous asymmetric vacuum. When this happens, 
the resulting ground state may have a crystal structure in charge density. (If the symmetric phase 
is unstable on the other hand, the resulting ground state would be a vortex lattice.) Then the 
translation symmetry must be spontaneously broken and there will be a sound wave as Goldstone 
boson. Fourth, the CTP symmetry is explicitly broken and so solitons and antisolitons have in 
general different mass spectrums. The Lorentz symmetry is also broken and some solitons can 
have zero or negative rest mass even though it must have positive kinetic mass. Fifth, the angular 
momentum operator from the Noether theorem should be modified. Some of these features are 
studied previously for different models [^, ^, 

In this paper, we concentrate on the self-dual Chern-Simons Higgs models with a uniform 
background charge density. As we shall see shortly, the usual extended supersymmetry associated 
with many self-dual models is no longer manifest in the presence of the background charge density. 

We start by introducing self-dual models in Sec. [2] , and then go on to study the vacuum structure 
of the model in Sec. [3]. We study the homogeneous symmetric vacuum that exists for all self-dual 
couplings, and find it stable even when the mass of elementary excitations are imaginary. The 
vacuum structure of the asymmetric phase depends on the parameters of the model and turns out 
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to be quite rich. We repeat the stabihty analysis for each homogeneous asymmetric vacuum. Rather 
interesting species of self-dual solitons appear in some of these vacua, which will be investigated 
in Sec. [4]. Here, we will find that the symmetric phase is in fact infinitely degenerate for certain 
parameter range. The section [5] is devoted to the matter of conserved angular momentum. The 
conserved Noether angular momentum turns out to be inappropriate due to severe divergences, 
and needs to be modified. A satisfactory and physically well-motivated finite expression is found 
by considering the same system on a sphere. Some concluding remarks are discussed in Sec. [6]. 

2 Model 

The first self-dual model, which is the one we will study in detail here, is the theory of a Higgs field 
(j) = fe^^ coupled to a Chern-Simons gauge field A^. The gauge field is coupled to a uniform 
background electric charge density so that the Lagrangian of the model can be written as, 

C='^e^-PA^,d,Ap-r\D^4>\^ -V - peA^, (1) 

where the self-dual potential is given by a specific form: 

u = ^\^\\m'-v'?-'^{m'-v\ (2) 

with D^(\) = {dfj_ + iA^)cl). The parameters v'^ and pe can be either positive or negative, but without 
loss of generality we may assume that k > 0. By shifting Ai ^ Ai — Pe^ijX^ /2k, we can see that 
the background electric charge pe is equivalent to a uniform background magnetic field = Pe/ f^- 
The C, P, and T are all broken: the Chern-Simons term breaks the P and T, while the background 
charge term breaks the C and CTP transformation. 

The second model is a Maxwell Chern-Simons Higgs theory with a neutral scalar field N , whose 
Lagrangian is 

^ = "4^^^'^^''' + ^e^^^M^S^ylp + -^{d^^Nf + \D^(p\'' -U- p^Ao, (3) 
where the potential is 

U = iV^l^p + ^{2\^\' -v^- 2KNf - peN. (4) 
o 

Actually, the first model can be obtained from the second by taking the limit ^ 00. The Maxwell 
term and the kinetic term for become negligible, while the field equation for N may be used to 
recover C from C. Similar models without the background charge density have been studied before 
[^, and also the pure Maxwell-Higgs system with k = was studied extensively by one of the 
authors |Q. Except for the case k = where the parity is a good symmetry, there is no qualitative 
difference between the two models and so we will focus on the simpler pure Chern-Simons-Higgs 
model C. 
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The system is invariant under the local gauge transformations. Gauss's constraint obtained by 
varying Aq is 

KFi2 + fie + Ao)-pe = 0, (5) 

where the dot denotes the time derivative. The electric charge current —f'^{d^6 + A^) is conserved 
and the corresponding conserved charge is Q = — J (Pxf'^{6 + Aq). 

Note that both the total electric charge and the total magnetic flux must be conserved inde- 
pendently, satisfying Gauss's constraint. (One can imagine the theory is on a spatially compact 
manifold so that "total" magnetic flux and "total" electric charge make sense.) Therefore, we 
find infinitely many superselection sectors, each of which is labeled by the total electric charge 
or magnetic flux of the system. Later on, we shall consider cases where the background charge 
density is cancelled at spatial infinity by either magnetic flux F12 or the Higgs charge density 
—f'^{9 + ^o)- Collectively, the first class of sectors are to be called the symmetric phase and the 
latter the asymmetric phase, respectively, for the obvious reason. 

Since the background charge density is uniform, the space-time translation symmetry is pre- 
served. The corresponding energy-momentum tensor T^,y obtained through the Noether theorem 
is 

T^. = ^e^paA'^d.AP + d^fd.f + fid^e + A^)d,9 - 7]^,C. (6) 

The conserved energy of the system is then E = JcPxTqq. After a partial integration. Gauss's 
constraint (||) can be used to show that 

E = l d'x[lp + + \ne + A,f + + A,f + [/}. (7) 

By using Gauss's constraint (|5|) again and integrating by parts, we rewrite the energy as 

E = j d^^[\f^ + + '^ifi9j0 + A,)f + ^-f{e + Ao-^{f- v')f^ 

+ yJ ^'^^12 - J d^x d^[^ei,f{dje + Aj)]. (8) 

The last term vanishes identically, as long as F12 = or = asymptotically. Recall that these 
conditions are characteristic of the asymmetric phase and the symmetric phase alluded to earlier. 

In the asymmetric phase, we naturally consider the excited states of finite total magnetic flux 
^ = J (PxFi2, and for those configurations we find the following self-dual bound, 

E > -M/. (9) 

As we will see later, within certain parameter range, this bound can be saturated even when 
is actually negative. 

In the symmetric phase, however, there exists a finite average magnetic flux density, F12 = Pe/ 1^- 
Then, the symmetric phase comes with a nonzero vacuum energy, Eq = {v^ /2) J F^2j is evident 
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in Eq. (|^. On top of this, of course, there will be locaHzed net fluxes of soUtonic configurations. We 
must then consider the excited states of finite net excess mae; netic flux ^' - ^'o = / d^x {Fu - Ffa), 
the excess energy of which is again bounded below 

2 2 

E-Eoy^i"^- ^o) = -y / d^^fio + Ao). (10) 

In order to saturate these self-dual bounds, the modulus field / must be static and, together 
with the gauge field A, solve the following set of the first order self-dual equations: 

dif + eijf{d,e + Aj) = (11) 

KF,2 + ^f\f-v^)-Pe = (12) 

Here we have used Gauss' law to remove + Aq in favor of k.Fi2- 

Are there any vortex-like configurations solving these self-dual equations? First of all, the 
modulus field / must vanish at the center of a vortex in order for the Higgs field to be well-defined, 
and this combined with Eq. ([Tl] ) implies that a self-dual vortex-like configuration always consists 
of anti- vortices only, 9 = — X^a^^S (^ ~ Qa) + V where ija are the positions of the anti- vortices and 
is a single valued function. 

Then, we may combine the coupled first-order equations above to produce a single second order 
equation for / with sources at the sites of anti- vortices. 

VHnf-^[f{f-v^)-2Kpe]=i7rY,6ix-qa). (13) 

a 

In the symmetric phase, for instance, we see that, far from the solitons, the modulus field behaves 
exponentially / oc q-Pb^'^/'^'^ so that there is no self-dual configuration in the symmetric phase 
whenever pe < 0- 

(The solutions of this last equation in case the external charge density vanishes were studied in 
detail before [^. In the symmetric phase, there are q-balls and q-balls with vortices, while, in the 
asymmetric phase, there are topological vortices. One of interesting properties of these solitons is 
that they carry the fractional spin and satisfy the fractional spin statistics therein.) 

3 Vacuum Structure 

Next, let us consider the vacuum structure. Of particular interest are those superselection sectors 
with homogeneous vacuum states in it. One may naively expect to find a homogeneous ground 
state in any given sector, but it is not difficult to see that there exist only two cases where this is 
possible: either J P {0 + Aq) = or J = 0. 

The homogeneous ground state one finds in the first sector is just the uniform magnetic field 
configuration, Ffg = Ps/k, without any Higgs field expectation value, (</>) = 0. Since this configu- 
ration saturates the self-dual equations trivially, this symmetric homogeneous vacuum must be an 
absolute minimum within this superselection sector. 
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Prom the form of the self-dual potential U{f), however, one can see that for sufficiently large 
and positive pe, the potential is concave at origin = //\/2 = 0, and the homogeneous symmetric 
vacuum appears to be at a saddle point, unstable even perturbatively. What are we missing? The 
uniform magnetic field Pe/i^ affects the dynamics of the small perturbation profoundly, among other 
things. To the first order in perturbation, the small deviation 6(j) = 4> obeys the following equation 
with a nontrivial kinetic part.. 

-dl5<l) = -{di - iA,if 5^ + U"{0) 5(^, (14) 

where the gauge field configuration A^i satisfies 51^4*2 — d2A^i = pe/n. To this leading order, the 
gauge field fiuctuation A — A^: decouples from the Higgs mode, and actually vanishes identically. 

Note that a perturbative instability is possible if and only if the operator on the right hand side 
of Eq. (14) possesses a negative eigenvalue. But, this eigenvalue problem is just the well-known 



Landau level problem with an extra "potental" term U"{0). The eigenvalue spectrum of the kinetic 
part is (2n -|- 1) |pe|/'« for all integers n > 0, which implies that the operator is bounded below by 
\pe\/i^ + U"{0) = \pe\/i^ + [v'^/Ak^ — Pe/n] > 0. Hencc, the perturbative calculation predicts that 
the homogeneous symmetric vacua is indeed linearly stable. By the way, the infinite number of zero 
modes that appear when v'^ = (with positive pe) is intimately related to the vanishing self-dual 
energy-bound @. 

In the second sector, where Gauss's constraint is satisfied by introducing an opposite scalar 
charge density, it is again possible to obtain a homogeneous configuration that solves the field 
equations. With F12 = everywhere, the constraint is then satisfied as a(/) = pe/f^ = + Aq. 
The eff'ective energy density of such configurations is then given by the sum of the scalar potential 
and an "electrostatic" contribution coming from the scalar kinetic term. 

£ = dl + U{f) = if - v') - 2Kp,f. (15) 



Searching for homogeneous vacua in this second super selection sector is now simply a matter of 
minimizing £ with respect to /. 

Once we find a local minimum at / = u, the next logical step is to test its stability. One 
convenient way is to check whether the configuration saturates the self-dual bound. Since there is 
no net magnetic fiux, the self-dual bound of the energy functional vanishes, and a minimum of £ 
saturates the bound if and only if the value of £ there is identically zero. Later on, we will find 
that this does not necessarily hold for all available minima. 

On the other hand, we also have the option of carrying out a perturbation to test the linear 
stability. For this, one again expands the field equations around the homogeneous vacuum. But 
unlike the previous case of symmetric vacuum, the gauge field fluctuation does not decouple, and 
we find a system of coupled linear partial differential equations with uniform coefficients: 

[-dodo + didi + a^{u)-U"{u)]5f + 2ua{u)5Ao = 0, 
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K eij di 6Aj + 6Aq + 2u a{u) 6f = 0, 
K eij {do 6Aj — dj dA^) + 6Ai = 0, 



(16) 



where we used the gauge = 0. A mode expansion with the frequency w and the spatial momentum 
p yields the following dispersion relation: 



1 



w 



+ 16a2(n) 



(17) 



Since £" is necessarily nonnegative at local minima, large wavelength fluctuations are massive ones 
with the masses squared given by 

,4 



u 

72' 



(18) 



These are associated with small fluctuations in the Higgs mode and in the massive gauge field 
mode respectively. When P(, = 0, the spin of the Higgs particle of mass itlh is zero, while that of 
the vector particle of mass rriA is one. By the continuity, we expect this aspect persists even with 
nonzero charge background. At the zero momentum, Eq. ( [T^ ) implies that the Higgs mode and 
the vector field mode are decoupled, confirming that the spin of the Higgs particles and the vector 
bosons are zero and one, respectively. We also note that two branches of the dispersion relation do 
not cross each other except at = when rriA = mriH- 

A linear instability occurs if and only if w'^ becomes negative at some nonnegative value of . 
We find the homogeneous vacuum of / = u is linearly unstable if and only if the following inequality 
is satisfied. 

Aa^>{mH + mAf. (19) 

This can be easily seen from the explicit forms of w"^ and when w'^ takes the minimum value, 
given the dispersion relation above: 

1 



16a2 
1 

16a2 



16a" 



{ruH - ruA 



4a' 



(20) 



The w'^ can take a negative value if either Aa'^ > (m// + m^)^ or {mn — itia)'^ > Aa?, but the latter 
implies pi < 0, and thus is unphysical. 



In some cases, there exists a so-called magnetoroton mode that are excitations with the least 
energy but at nonzero momentum. In the present context, magnetorotons should appear whnever 
p1 > 0. Depending on whether mfj > niA or not, the roton mode occurs on the branch of the Higgs 
mode or the vector boson mode. One thus naively expect that the spin of the roton mode is zero or 
one, depending on whether mn > tua or not. However, there is no rest frame for rotons and it is 
not clear at this moment whether it is possible to separate the orbital and spin angular momentum 
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for rotons. As we shall see shortly, magnetorotons in a stable vacuum exist for a limited range of 
Kpe and that only for positive v^. 

Now let us list all possible minima of £ over various ranges of the parameters v'^ and np^- 
Extremizing £, we find the following four solutions: 

u± = ^ , 

4 ^ -'^V-'-^^-Pe, (21) 

Not all of these represent a physical vacuum. Depending on the parameters of the theory, some of 
them become negative or even complex. Also, a > may correspond to a local maximum instead 
of a minimum. Through a careful but elementary study of £, one finds the following homogeneous 
vacuum structure: 



(1) < 0, < npe] u-i- is an absolute minimum and saturates the self-dual bound. 

(2) t>2 < 0, Kpe < 0; n_|_ is an absolute minimum but does not saturate the self-dual bound. 

(3) v'^ > 0, npe > r;^/24; is an absolute minimum and saturates the self-dual bound. 

(4) > 0, v^/24: > upe > 0; u-i- is an absolute minimum and saturates the self-dual bound; n_ is 
a local minimum that does not saturate the self-dual bound. 

(5) > 0, > Kpf, > —v'^/8; u± are two degenerate absolute minima and saturates the self-dual 
bound. 

(6) > 0, — > upe] u+ is an absolute minimum but does not saturate the self-dual bound. 



The cases where the equality rather than the inequality holds can be understood as the degenerate 
limits of the cases in the above list. 

Despite this long list of different cases, the linear stabilities are more or less determined by the 
vacuum expectation value u alone. For example, the homogeneous vacua u = u±, whenever they 
are real, do saturate the self-dual bound £ = 0, implying that they are always stable under linear 
perturbations. Indeed, an explicit calculation shows that 

4a^(n+) = {mH{u+) - mA{u+))'^, 

4a^(n_) = {mH{u^) + mAiu+))^. (22) 

Clearly both u± vacua are linearly stable, in view of the instability criterion ([l9|) . The above 
equation implies p^(n+) < so that there is no roton in the u+ vacuum. Interestingly enough, u- 
vacuum is actually right at the borderline of instability, and generically possesses a massless roton 
mode at p1{u-) = niHiu-) mA{u-). 

The linear fluctuations around u±- vacua are a bit more complicated. First of all, neither 
nor n_ saturates the self-dual bound since £{u±) > unless Kpe = or Kpe = — f^/8: Even the 
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linear instability becomes a nontrivial issue. Using Eq. (21), we find the following behaviour of wl 
and pI in terms of the vacuum expectation values u±: 



k2 {v^ - 3n^)2 



2 _ {2ul -v^) 

where we have traded off pe in favor of the vacuum expectation value u±. 

In u+ vacua for both cases (2) and (6), it is easy to see that w'^ is always nonnegative, which 
again translates into the linear stability. Are there roton modes? The answer is yes, but only for 
case (6). pi is positive in this case, provided that -u^ < f^. In terms of Kp^,, the corresponding 
range is given by > Kpe > —v^. 

Combined with the case (5) where a roton mode is found for n_ vacuum, this means that there 
exist a roton mode around a homogeneous asymmetric vacuum only for theories with > and 
> Kpe > — v^. This roton is massless in u- vacuum while massive in n_|- vacuum. Note that the 
roton mode can occur on either branches of the dispersion relation. 

The only remaining vacuum to consider is that of the local minimum = that exists for 

> and > Kp^ > 0, namely for the case (4). Although it is a local minimum of the effective 

energy density £, it is always unstable. Indeed, as we vary upe within this range, interpolates 
between and 0, and Eq. ( p3|) then implies that w1 is always negative while pi remains positive. 
Hence, this particular local minimum is classically unstable against an inhomogeneous fluctuation 
(of length scale ~ k/vu-), unlike any other case we considered before. At this point it is not clear 
whether such instability will lead to excitations on the homogeneous self-dual vacuum = u^. 

4 Self-Dual Solitons 

Having studied the homogeneous vacuum structure in great detail, we are now in position to ask 
what are possible self-dual soliton solutions in the symmetric and the asymmetric phases. Previous 
works shows that in the models without the background charge density there can be q-balls in the 
symmetric phase and vortices in the asymmetric phase Q. Do they persist in our case? If so, 
how does the background deform them? Are there other kind of solitons? Specifically, we must 
investigate possible self-dual configurations that asymptotically approach a given homogeneous 
vacuum. 

Similarly to what we have seen in the study of homogeneous asymmetric vacua, the answer 
depends on the parameters of the theory and also on which homogeneous vacuum we consider. For 



instance, in u± vacua, the second self-dual equation ([1^) implies a divergent total magnetic flux, 
thus making a self-dual soliton impossible. Another case where a self-dual soliton does not exist is 
the symmetric phase in theories with negative pe, as pointed out earlier. 
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Are there self-dual solitons in u± vacua? While a generic self-dual configuration would not 
be rotationally symmetric, let us focus on the rotationally symmetric solutions to gain an insight 
on the qualitative features of the self-dual solutions. The ansatz for a rotationally symmetric 
configuration around the origin is made of /(r), 6 = —mp, Ai = A{r)(pi. Rewriting the second order 
form of the self-dual equation in terms of a new field y = log(/^/|t'^|) and the dimensionless radial 
coordinate s = r|f^|/«;, we find a second order dissipative equation for rotationally symmetric 
self-dual configurations: 

ds'^ s ds dy s ' 

where y = ln(/^/|t;^|) and V{y) = — ±6^ + 2Kp(,y/v'^, ±1 being the sign of v'^. 

Let us first consider the u+ vacuum. Note that y^ = ln(n^/|f^|) is always a local maximum 
of V{y)- In Eq. (24), a self-dual anti-vortex in the u+ vacuum starts at the asymptotic value 
y-^. = log(ii^/|t'^|) and smoothly evolves to y = —oo, as r decreases from oo to 0. The dissipative 
term tends to increase the total "energy" as r approaches the origin, and thus lets = |u^|e^ 
reach zero at the center of the vortex even when V has a linearly increasing potential barrier 
toward y = —oo {pe < 0). 

However, if the system starts out at the other possible asymptotic value y_ = log(n?^/|f^|) that 
exists in the case (5), the only nontrivial solutions at large distances exhibit the following damped 
oscillatory behaviour: 

oc ^e^''^^"^y-l (25) 
Vs 

Although this lets settle down to at large spatial distances, the resulting configuration does 
not have a convergent total magnetic flux, 

r dr rFi2 ~ y/^e^"'^'"^^/"" (26) 
Jo 

If we are considering a rotationally nonsymmetric soliton instead, that will simply introduce an 
extra term d'^f / in the equation above that cannot change this asymptotic behaviour. There is 
simply no nontrivial self-dual configuration that asymptotes to the homogeneous vacuum u_. This 
is clearly related to the fact that there are massless roton modes in this vacuum. 

This leaves us with only two choices for the asymptotic state of self-dual solitons: either the 
homogeneous asymmetric vacuum = that exists when Supf, -|- v'^ > or the the homogeneous 
symmetric vacuum = with Pe ^ 0. The second half of this section is mainly devoted to 
understanding novel solitonic states that arise in these vacua when the background pe is nontrivial. 
As we will see shortly, not all of the self-dual solitons here survive the limit pe 0. 

First of all, in the asymmetric phase with v'^ > 0, the self-dual solitons are rather similar to 
the zero background case. Figure 1. shows the self-dual anti- vortex configuration for n = 1, > 
and Kpf. = v"^. 
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Figure 1 : Plots of a typical rotationally symmetric anti- vortex in the asymmetric phase that 
saturates the self-dual bound. The first graph shows (solid line) and k.FyiI Pe (broken 

line) as functions of the distance from the center. The second graph is a plot of the flux (divided 
by 27r) within finite distances from the center. The radial distance is measured in terms of k/ | | . 

These self-dual anti-vortices do survive the limit pe 0, and become the ordinary topological 
solitons of the background-free theory. However, if happens to be negative, an unexpected new 
type of solitons arises, as we discuss below. 

One interesting aspect of the self-dual bound here is that the energy bound is proportional 
to the magnetic flux rather than to the absolute value thereof. As a result, a unit self-dual anti- 
vortex in the u+ vacuum of case (1) actually has a negative energy — 7r|u^|. Of course, this does 
not mean that the homogeneous vacuum of case (1) is unstable against the proliferation of anti- 
vortices, for whenever an anti-vortex is created, a vortex must be also created to conserve the 
total magnetic flux. In fact, as mentioned before Eq. (p!3|), and due to the broken OTP, a vortex 
does not saturate its self-dual bound 7r|t;^| > 0, for it cannot solve the self-dual equations at the 
vortex center. Hence, a vortex-anti-vortex pair always costs more energy than 7r|w^| — ■t:\iP'\ = 0, 
and the homogeneous asymmetric vacuum is stable against this particular process. Even though 
anti-vortices have negative rest energy, we see their kinetic mass must be positive since the self-dual 
energy bound implies that slowly moving anti-vortices should have a positive kinetic energy. 

Still, one may wonder what happens to such negative energy solitons as we take the limit /?e — > 
where the self-dual bound of the soliton energy is clearly positive. A useful equation to look at in 
order to understand what happens, is the second self-dual equation ([l2[) that can be solved for 
algebraically: 




(27) 
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The flux density F12 must be bounded above by pe/i^ to maintain real /, while the total flux of the 
anti-vortex is quantized at a positive value / F12 = 2tt. Therefore, as pe/i^ — > 0, the soliton core 
must be spread out to infinity, leaving behind a symmetric vacuum, = 0. 

What about the symmetric phase? Are there any negative energy q-balls, possibly with an 
anti-vortex embedded inside. Again using the the second self-dual equation, we can express the 
energy bound as follows, 

E-Eo>^-^{^-^o) = "-^Jd'x{Fu-'-^) = ^,Jd'xf{v'-f). (28) 

For negative v'^, the right-hand-side is manifestly positive, while the same could be true for positive 
v'^ provided that remains smaller than v'^. 

However, it is easy to see that, for a strictly positive Pf>, may grow larger than f^. To 
see this, again consider Eq. (p^), where the "potential" energy V now has a unique maximum 
at ?/+ = log(u^/w^) > but no minimum. In this picture, a self-dual q-ball would interpolate 
smoothly between y = —00 at r = 00 and some yo at r = 0, without the delta function 

source on the right-hand side. The only restriction on yo is that it be on the left hand side of 
the maximum: yo < which translates into /^(O) < u^. Then, a q-ball may have a central 
region (u^ < < u\) of negative charge density (F12 — Pe/^ < 0), surrounded by an outer region 
(0 < < v'^) that has a positive charge cloud {F12 — Pe/n ^ 0). In figure 2, two typical shapes of 
q-balls are depicted, depending on whether /^(O) is smaller or larger than v'^. In all figures in this 
section, the parameters are chosen such that Kp^/v^ = 1. 




Figure 2: Plots of two q-balls in the symmetric phase with npe/v'^ = 1. Solid lines are /^/w^, 
and the broken lines are the excess magnetic field, or equivalently, the Higgs charge density 
divided by k, in the unit of v'^/k^. The first graph corresponds to /^(O) = 0.9 while the 
second to p{0) = lAv'^. 
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Furthermore, as /^(O) u\, the central region of negative charge density becomes larger and 
larger, so that, in the strict limit /^(O) = u\, it overtakes the entire plane such that the correspond- 
ing self-dual configuration is simply that of the homogeneous asymmetric vacuum f^{r) = u\ with 
Fi2 = 0. In other words, there exist self-dual q-balls of arbitrary negative charge k — ^'o) < 0, 
or equivalently of arbitrary negative energy E — Eq = [v^ /2){^ — ^q) < 0. This is illustrated in 
figure 3. 




Figure 3: Plots of P/v^ and KFi2/Pe as the value of /^(O) approaches = 2u^. As /^(0)/u^ 
increases toward one from below, the central region of vanishing magnetic field grows steadily 
(the second graph). 

In the other limit of yo — oo, the charge cloud of the soliton remains manifestly positive 
everywhere, but the strength thereof decreases indefinitely so that the net flux also decreases 
indefinitely. The upshot is that we have a one-parameter family of rotationally symmetric q-balls 
such that their total charge is bounded from above by a finite positive quantity. Furthermore, for 
each net positive charge — ^'o) > allowed, there exist two different solitons, depending on 
whether /^(O) is smaller than or larger than v^. 

One interesting consequence of this is the degeneracy of the symmetric vacuum. For a moderate 
value of /^(0)/f^ > 1, there must exist a nontrivial configuration of zero total charge, thus of zero 
net magnetic flux and zero net energy. This is simply because there exists a one-parameter family 
of solutions that interpolate between the positive charge solution of (0) = v'^ and arbitrarily large 
negative charge solutions that show up as /^(O) u^. The zero-charge configuration looks like a 
small compact island of broken phase enclosed by a cloud of positive charge that exactly cancels 
the negative charge inside. In figure 4, such a configuration is depicted when npe/v* = 1. 
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Figure 4: Plots of a zero-energy solution. The second plot shows the net flux (divided by 2tt) 
within finite distances from the center of the configuration as a function of rv^/n. 



It is not too difficTilt to convince ourselves that similar considerations may be applied to q-balls 
with anti-vortices at the center. After fine-tuning the behaviour of near the origin for each fixed 
anti-vorticity n, we find that there exist rotationally symmetric solutions of arbitrary negative 
net energy E — Eq < 0, and a unique solution of zero net energy for each n. Such zero-energy 
configurations consist of a thick ring of negative charge sandwiched between a positive charge core 
and a surrounding outer region of positive charge cloud, such that the total charge is exactly zero. 

Hence, whenever both v'^ and pe are strictly positive, the homogeneous symmetric vacuum is 
degenerate with infinite number of rotationally symmetric self-dual configurations. Are there other, 
more complicated, states of vanishing net charge that are also degenerate with the homogeneous 
symmetric vacuum? One possible suggestion is that these zero-energy "solitons" do not interact 
with each other so that we may create many of them at once without costing any energy. This 
possibility along with many interesting questions concerning general multi-soliton configurations 
are postponed to a future study. 

5 Angular Momentum 

One of the more important aspects of Chern-Simon theories is the fractional spin and statistics 
carried by the solitons. In view of the fact that the background charge density alters the behaviour 
of the theory rather significantly, and in particular of the explicitly broken OTP, it would be very 
interesting to find out whether spin and statistics are modified. As a first step along this direction, 
let us consider the conserved angular (and linear) momenta. 
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The naive angular momentum density obtained by the Noether theorem contains a term, pro- 
portional to peAj, that not only cause divergences in all phases but is gauge-invariant only up to 
a total derivative. Explicitly, the Noether angular momentum is 

JNoether = " / eijx'l^fdjf + f{e + Ao){dje + Aj) - p,Aj^. (29) 

A conventional way of dealing with such problems is to isolate a finite part by adding an appropriate 
manifestly conserved current, and then justify that the final expression is a physically sensible one 
[Q. However, we want to present an alternate and certainly more elegant procedure by considering 
a similar system on a two-sphere of radius R in the limit R ^ oo. 

A two-sphere has isometrics in the form of 50(3) rotations, and there exist the corresponding 
conserved angular momentum 3- vector L. Although the naive Noether expression for the density 
contains a similarly problematic term as above, we may now perform a partial integration to recover 
a manifestly gauge-invariant expression for the density, since there is no longer a boundary to worry 
about: 

^ = - J R^dnl^fir X V)f + f{e + Ao)rx {VO + A) - p^r {r ■ V x A)|. (30) 

This satisfies the usual so(3) commutation relations. Here, r denotes the coordinate vector of the 
three dimensional Euclidean space where the spatial sphere is imbedded. In this notation, the 
actual spatial manifold corresponds to |r| = R. The 3-vector A must be tangent to this two- 
sphere, since it represents the spatial gauge field on the two-sphere. Finally, the total magnetic 
flux ^' = / R?dQ. r • V X A should be quantized in the unit of 27r to ensure the gauge invariance. 

Already, we can see the advantage of working on the sphere instead of the plane: the last term 
is now proportional to V x A, thus manifestly gauge- invariant. Furthermore, given a uniform 
radial magnetic field strength B, there is no BR'^ divergence in the total angular momentum L, 
due to cancellations between two hemispheres, unlike JNoether that has such a nasty divergence in 
the symmetric phase. In order to take the plane limit R oo, we imagine all interesting activities 
of the system occurs within finite distances from the north pole. 

First, let us consider the asymmetric phase where we have a uniform Higgs charge density 
—f'^{9 + Aq) = —pe at large spatial distances from the north pole. The vortices and the anti- 
vortices of total flux ^' = 2re7r, if distributed near the north pole, cause a quadratic divergence 
^ PeR^ in the third component L^. However, this should not be surprising at all. A similar 
divergence arises if we consider an electrically charged particle on a large two-sphere threaded by 
a uniform magnetic field B of fixed strength. The extra angular momentum ~ egr, familar from 
the studies of charge-monopole system would diverges quadratically simply because the total 
magnetic charge g = BR^ diverges quadratically. In the asymmetric phase, there exists the so- 
called Magnus force, between the uniform Higgs charge density —pe and the localized fluxes, which 
simulates the Lorentz force. 

As long as we are concerned with the plane limit, we may as well remove this divergence in each 
superselection sector of fixed ^. Then, we find the following modified angular momentum, valid 
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for the symmetric phase on the plane: 

fd,f + f\e + Ao)idje + Aj] 



+ ^^'F,,\. (31) 



Note that the vestige of the Magnus force contribution still remains in the last term. To see this 
more clearly, consider a generic 9 field configuration, 

9 = -J2{-l)aATg{r-qa{t))+rj, (32) 

a 

where rj is, as usual, single- valued. Then, the above angular momentum can be rewritten, after an 
integration by part, as 

Jasym = - f(fx X^Eij {fd,f + [f {9 + ^o) " Pe]{dj9 + Aj)} - TTp,J2{-l)a\qa\^ . (33) 

The last term that sums over the locations of vortices, is similar to what one would get for a charged 
particle moving under a uniform magnetic field. Vortices in our case feel the aforementioned Magnus 
force instead, and this term represents the corresponding modification. For a similar expression in 
the pure Maxwell-Higgs systems, see Ref.|||] 

In the large sphere limit, the linear momenta on the plane are found from its relation to the 
SO (3) angular momentum via = —RP^ and = RP^. Note that as long as the nontrivial 
configurations are concentrated near the north pole there will be no quadratic divergences in and 
L^. Alternately, we may extract P*'s from the plane angular momentum, for the latter transforms 
like J — > J + eijtt^P^ under an infinitesimal translation x* — > + a*. The two resulting expressions 
coincide with each other. Prom the quadratically divergent part of the commutation relation 
[L^,L'^] = iLp' in the plane limit, we find [P^,P^] = ipe'^- This modified form of the commutator 
is again due to the Magnus force exerted by the Higgs charge density —pe- Here we have to use the 
commutation relation [Ai{x,t), A2{y.,t)\ = i5'^{x — y)/^. 

In the symmetric phase, the field / will be nonzero only near the north pole so that at large 
distances from the north pole, there is a uniform magnetic field f • V x A = pe/^- Physically 
interesting configurations are then those with finite net magnetic flux, ^' — ^'o < oo where \E'o = 
AttR? Pf./ n. Again nontrivial configurations are assumed to be within finite distance from the north 
pole as R ^ oo. In order to take the plane limit, it is sensible to subtract = R^ J dO, P(,r, and make 
the angular momentum density to be concentrated near the north pole. This effectively replaces 
f-VxAbyr-VxA — pe/K in the last term in Eq. (|30|). 

Now, diverges quadratically in the plane limit as R^pe (^' — ^o)- In the symmetric phase, 
the origin of this divergence is even more transparent, for the system is basically that of charged 
particles (q-balls) of total charge k (^ — ^q); moving in a large two-sphere that surrounds total 
magnetic charge R'^pe/n at the center. In each superselection sector, we may again remove the 
quadratic divergences to obtain the modified angular momentum, valid for the symmetric phase on 
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the plane. 



Jsym = - j M/ + /'(^ + ^0)(9,^ + Aj) + (fiFu " /9e)}. (34) 



We can find the hnear momenta similarly as in the asymmetric phase, and can show that they satisfy 
a nontrivial commutator [P^,P^] = ipe{'^ — ^o)- This last expression may be rewritten, using 
Gauss's constraint, in terms of the Higgs charge density and the background magnetic field Pe/ n 
of the symmetric phase: [P-*^, P^] = —i {pe/ n) J d?x {6 + Aq). Written this way, the nontrivial 
commutator can be naturally attributed to the velocity-dependent Lorentz force on the q-balls, 
exerted by the homogeneous magnetic field pe/ k. 

Now that we derived the correct angular momenta in the presence of the uniform background 
fields, it is time to consider how these abstract expressions translates into angular momenta of 
specific solitons. For simplicity, we shall consider single static solitons. 

Let us start with topological anti-vortex in the asymmetric phase. Since we can safely exclude 
the center of the soliton at = from the space integration of Eq. (^) without affecting the 
value of the angular momentum, we may introduce a vector potential Aj = djO + Aj, and rewrite 
the Jasym as an integral over E^ = E? — {q]: 

Jasym = ~ j ^"^^ ^ij^' j/^j/ ~ f^AjekidkAi^ - n'npe\q\^. (35) 

For the static self-dual solutions (and also for rotationally symmetric solutions), diAi = and he 
integrand becomes a total divergence: 

Jasym = ^ j ^3 \^^^ ^k) - Aj{x'' Ak)^ - m:pe\q\^. (36) 

Note that ^, as a gauge invariant quantity, must be exponentially small at large distances so that 
the only boundary to speak of is at the center. If the anti- vortex has n topological charge (i.e., if 
the total magnetic flux is 27rn > 0), Aj near the center is dominated by djO ~ nejkX^ /r"^ , and the 
angular momentum is easily obtained as follows, 

Jasym = -VTKn^ + TTU {-pe) \q\^ . (37) 

In the limit pe — > 0, this reproduces the expected results found in Ref.|§, ^. The last term 
obviously reflects the Magnus force on the anti-vortex exerted by the background Higgs charge 
density — [9 + Aq) = —pe- Clearly, the above expression is true also for rotationally symmetric 
vortices. Thus, we have shown vortices and anti-vortices in Chern-Simons systems with nonzero 
magnetic flux carry nonzero spins — vr/tn^ that are independent of pe- 

In the self-dual case, we can have static multi-anti-vortex conflgurations where anti-vortices do 
not overlap each other. Eq. (5^) can then be explored much further to gain an understanding of 
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slowly moving vortices and thus of their statistics. (See Ref.[^ for the case pe = 0.) We hope to 
return to this topic in near future. 

In the symmetric phase, on the other hand, we introduce Ai = diO + Ai + eijx^ pe/{2K.). Using 
Gauss's constraint, we again find that the angular momentum for a static self-dual configurations 
can be written as an integral of a total divergence, 

Jsym = <freijx'Aj€MdkAi = d^x dj - . (38) 

As we have seen earlier in a more general context, the angular momentum in the symmetric phase 
must also contain a term ~ \q\'^ that originates from the Lorentz force exerted by the uniform 
magnetic field Ps/k.. However, this expression, compared to the analogue in the asymmetric phase, 
does not seem to contain such a term. The crucial observation that resolves this apparent puzzle is 
that, unlike A above in the asymmetric phase, A exhibits a powerlike behaviour at large distances 
so that we need to worry about the contribution from the boundary at infinity. In fact, if the 
soliton does not contain any vortex or anti-vortex at the center and thus is a pure q-ball, the only 
boundary is at infinity. 
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J,ym = i^f dl\-{A''Ak)-^—^\. (39) 



For a single, rotationally symmetric q-ball of net magnetic fiux lira, in particular, the self-dual 
equations (pi]) and (13) imply the following asymptotic behaviour. 



provided that the soliton is centered at = g*. After a careful usage of this limiting form, we find 
the total angular momentum of a self-dual q-ball of net magnetic fiux 27rQ at = g*: 

<^ym = 7rK|ap — vr/ta — (41) 

K 

The last term is exactly the position dependent term we need, which again shows that our angular 
momentum is sensible. 

In this section, we derived appropriate expressions of angular momentum in each of the sym- 
metric and asymmetric phases, and evaluated them on self-dual q-balls and anti- vortices. While the 
angular momentum acquires a new position dependent terms that refiects the interaction between 
the solitons and uniform fields of the ground states, it also include the usual fractional spin that 
is independent of the background. Those additional position dependent terms, which are propor- 
tional to Pe, arise because the slowly moving solitons feel a velocity-dependent force exerted by the 
uniform fields of the surrounding vacuum. 

6 Conclusion 

We studied the characteristics of a self-dual Chern-Simons Higgs model coupled to a background 
electric charge density. We found a rich vacuum structure, and subsequently tested the classical 
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stability of homogeneous vacua along magnetoroton modes. In certain stable vacua, novel self-dual 
solitonic configurations are found and studied in detail. Finally, the divergent Noether angular 
momentum is successfully modified to a sensible and finite expression. From this formula, we find 
that the soliton angular momentum is corrected by an effective Lorentz force, and acquires a term 
quadratic in the position as well as the usual (background-independent) fractional spins. 

As the structure of the theory is very rich, we were not able to cover even the classical aspects 
of the theory completely. For instance, there are cases where no self-dual solution exists: Although 
we still expect to find non-self-dual solitons in such cases, practically nothing is known about 
them other than the spin. Another important topic which is not covered here is the matter of 
the spin-statistics. Since our solitons carry nonzero fractional spin, we expect there is nontrivial 
statistical interaction between solitons. It would be interesting to find out whether the statistical 
interaction follows the naive expectation that it is made of the spin contribution and the background 
Magnus/Lorentz force contribution. The study of the classical dynamics of slowly moving anti- 
vortices in the asymmetric vacuum may lead to a better understanding in this regard. Finally, 
it is not clear to us at all how our self-dual model could be extended to have A'^ = 2 or = 3 
super symmetry, which is expected from self-dual models in general fTo|] . 

Quantum aspect of our theory should be also quite rich. Especially, the understanding of the 
infinite vacuum degeneracy of the symmetric phase poses a challenge. 
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